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Abstract 

We present a new, model-independent method to analyze radiative decays of mesons 
to a vector, isovector pair of pions of invariant mass square below the first significant 
7T7T threshold in the vector channel. It is based on a combination of chiral perturba- 
■ tion theory and dispersion theory. This allows for a controlled inclusion of resonance 

physics without the necessity to involve vector meson dominance explicitly. As an 
example, the method is applied to an analysis of the reactions n — > 7r + 7r~7 and 



7/ —7- 7T + 7T 7. 



1. Radiative decays are known to be very sensitive tools to explore decay 
mechanisms. Especially, when studied together with two hadrons as decay 
products, they enable us to adjust the invariant mass of the two-hadron system 
via a variation of the photon energy without interference of strong three-body 
final state interactions. 

We here present a new, model-independent method to analyze the mentioned 
radiative decays. The method is based on a combination of chiral perturbation 
theory (ChPT) - extended from SU(3) to U(3) [T|2j - and a dispersive anal- 
ysis. It is general and can be applied to all decays of mesons with 7T7T7 final 
states, where soft bremsstrahlung does not occur and where the pion pair is 
of invariant mass square below the first significant im threshold, which is, in 
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the (iso)vector case, the ujtt threshold. In this work, however, we focus on the 
decay of a pseudoscalar, i.e. the 77 and rj' meson, to a photon and a charged 
pion pair. In this case the general selection rules enforce the pion pair to be in 
the isovector channel. Especially for the decay of a relatively light meson, like 
the 77, one may expect that the ideal analysis tool is the effective field theory 
for the standard model at low energies, ChPT, and indeed the corresponding 
calculation has been available to one-loop order for a long time [3] . But com- 
pared to modern data one observes a significant deviation between the theory 
predictions and data — the source of which is mostly the non-perturbative 
mi final state interaction. Several efforts have been made to include final state 
interactions by unitarized extensions to the box-anomaly term - the latter 
of course determines the 77 — > txtx^i decay in the chiral limit. The tree- level 
calculation can be enhanced by a momentum dependent Vector Meson Domi- 
nance model [1] or more elaborate calculations in the context of Hidden Local 
Symmetries [5] . On top of the results at the one-loop level an Omnes-function 
can be applied to describe the effects of p-wave pion scattering [6fT] (see also 
[8]), or as done in the Chiral Unitary approach, a Bethe-Salpeter equation 
with coupled channels can be used to generate resonances dynamically [9] . We 
will later specify the details which discriminate these mechanisms from our 
method. The general problem with the vector meson dominance model is a 
priori that it is unclear what relative strength is to be put between the tree 
level contribution and the resonance contribution. This problem is resolved in 
the dispersion theory approach as we will see below. 

For the transitions at hand it appears necessary to disentangle perturbative 
and non-perturbative effects in a controlled way. The method is therefore split 
in two steps. In the first step the spectral decay data are fitted with a function 
of the form (the details will be given below) 



where the normalization parameter A has the dimension of mass 3 and where 



collects phase-space terms and the kinematics of the absolute square of the 
simplest gauge invariant matrix element (for point-particles). The latter is 



expressed through the 7T7r-two-body phase-space ai^s^) = J 1 — Ami/ in 
terms of the invariant mass square s n7T of the pion pair, while m (m^) denotes 
the mass of the decaying particle (charged pion). Since the initial state is a 
pseudoscalar and the final state contains a photon, the partial wave of the 
charged pion pair is expected to be dominated by 1 . If it can be confirmed 
by experiment that the other partial waves can be neglected, then the factor- 
ization given in Eq. ([T]) is exact and can be straight forwardly derived using 
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dispersion theory (see next section). If, however, it eventually might turn out 
that higher partial waves would be needed to get a precise fit of the angular 
data, then their contributions could still be added in a perturbative way to 
our model-independent expression for the p-wave decay amplitude. 

The pion vector form factor Fv(s nn ) is known very well from both direct 
measurements of e + e~ —> [Toll 11 11121113] as well as theoretical stud- 

ies |14|15II16II17II18II19II20] . It collects all non-perturbative irir interactions and 
is universal. On the other hand, the function P(s 7T7T ) as well as the normaliza- 
tion factor A are reaction specific and - at least for the decay of light mesons 
or, more accurately, for small values of - are expected to be perturbative 
in the sense of ChPT. In case of the 77 and 7/ decays in the focus here, left-hand 
cut contributions are suppressed both kinematically, since the particle pairs in 
the t-channel are to be (at least) in a p-wave, and chirally, since the p-wave 
7rr) interaction starts at next-to-leading order only [2~TH2"2"] . We may therefore 
expand P(s nn ) in a Taylor series around s nn = and define 

P(s^) = l + as^ + 0(sl 7 ) . (2) 

At higher order non-analytic terms, mainly from left-hand cuts, need to be 
considered. The parameters a and A allow insights into the physics underlying 
the decay process. Thus, in the second step of our method a proper matching 
scheme is formulated, to relate the parameters A and a to the parameters 
of the underlying effective field theory. For the example at hand we will find 
that this matching allows us, under certain assumptions, to impose a relation 
between t] — > 7r + 7r~7 and r]' — > 7r + 7r~7. 

The paper is structured as follows: after a brief discussion of the pion vector 
form factor, we apply the formalism to the mentioned r\ and 7/ decays in Sec. 3, 
extracting the phenomenological parameters and A^ . In Sec. 4 those are 
then interpreted via a matching to one-loop U(3) extended ChPT. In Sec. 5 an 
interpretation of essentially half of the empirical value of is given, whereas 
Sec. 6 contains a comparison to earlier studies. We close with a summary. 

2. We start with a brief discussion of the pion vector form factor. In terms of 

the vector-isovector current VH, it is defined via 

fj, 

(Tr + (p f )n-(p)\V*\0) = (p-pXF v (s W7T ) . (3) 

In the elastic regime the form factor is, in symbolic notation, defined via 

F v = M v + T^G^My , (4) 

with My, G n7T and T nn for the production vertex, the two-pion propagator 
and the tttt scattering amplitude, respectively. For simplicity, in this work we 
assume My to be real which is exact for the transitions we will study explicitly 
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below. Vector meson dominance models typically model the form factor, either, 
in its simplest variant by a single term, Fy(s nn ) = —m 2 p / (s n7T — m 2 + im p Y p ) 
where m p and T p are the p mass and width, respectively, or by writing the 
second term of Eq. (jlj) as a (sum of) vector meson propagator (s) times s n7T , 
which leaves the relative strength of the first and the second term as free 
parameter. We will here take a different route which leaves no freedom of 
choice. From the definition of Eq. (jlj) it follows straight forwardly that 

lm(F v (s n7T )) = <x( 

This relation holds for the whole elastic regime which, in case of the pion vector 
form factor, extends to values of s nn well beyond 1 GeV 2 , although already at 
s n7T = 16m 2 formally the first inelasticity opens. Eq. (jSJ is one way to present 
the Watson theorem (23]: since lm(Fv) is a real quantity, the phase of the form 
factor has to agree with the phase of the elastic scattering amplitude. In the 
elastic regime this equation can be written as lm(F) = tan<5 n Re(-F), with <5 n 
for the elastic tttt phase shift in the vector channel. The resulting dispersion 
integral can be solved analytically to give 



F v (s nn ) = exp (\s^(r 2 ) + ^ f°° ds _. 5ll{s) , . ) . (6) 

V TT Jiml S 2 (S - S w7r - 16) j 

Contrary to the standard procedure [IB] , we follow Refs. [H?f2"0"] and use a twice 
subtracted dispersion integral in order to guarantee that the integral over the 
phase converges in the elastic regime. The phase 5n can be taken from data 
or from theoretical analysis. The phase used in the present analysis is taken 
from Ref. [T7] (see Eqs. (A7) and (A8) therein), valid up to y/s cu t = 1.4 GeV 

- we use a smooth extrapolation of the phases beyond this energy following 
Ref. [H] . The phases agree with those of Ref. [21] up to 800 MeV and with the 
available data. Once the phase is fixed, Eq. fl6]) has only one free parameter 

- the subtraction constant (r 2 ), the mean square charge radius of the pion. 
Apart from the region around s n7T = m^, where p — u mixing shows up very 
prominently, we find an excellent fit of the data for (r 2 ) = 0.437(3) fm 2 . The 
uncertainty includes the uncertainty in the tttt phase shifts given in Ref. [T7] 
as well as the weak dependence on s cut . The form factor parametrization, 
including the uncertainty, is shown as the red band in Fig.[TJ The range for 
the radius is consistent with (r 2 ) = 0.452(13) fm 2 of Ref. [25] and (r 2 ) = 
0.435(2) fm 2 of Ref. [IB]. The possible impact of p — u mixing on the t]' decay 
spectra is briefly discussed below. 

As it provides the explicit solution for Fy defined in Eq. (T4]) , Eq. (E]) contain 
both the Born term (pions going out without interaction) as well as the final 
state interaction. In general, the mentioned dispersion integral fixes the form 
factor up to a multiplicative function that does not have right-hand cuts. 
Since the right-hand discontinuity of the transition amplitude t/') — > 7r + 7r - 7 
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Fig. 1. The (red) solid band shows the form factor derived from Eq. ([6|), the (blue) 
dashed line the result from one-loop ChPT (see Eq. ([7])) — both with identical values 
for the pion radius. The time-like data from Refs. [11] and [12] are shown as solid and 
open circles, respectively. The space-like data are from Ref. [13]. The short (long) 
thick, horizontal bar in the left panel denotes the kinematic range covered in the 
decay of the 77 (77'). The right panel shows, as a linear plot, a zoom into the s n7T range 
relevant for the r/ decay. In this energy range the form factor can be approximated by 
the polynomial \F v (s n7T )\ » l + (2.12±0.01)* w „. + (2.13±0.01)a^ r + (13.80±0.14)sj w 
with s n7r measured in units of GeV 2 . 

agrees with that of the pion vector form factor, the factorization employed 
in Eq. ([T|) is justified — under the above-stated qualifications that the higher 
partial waves can be neglected and that the left-hand cut contributions are 
suppressed. We chose the standard normalization for the form factor, namely 
Fy(0)=l, which corresponds to My = 1 in Eq. (0|. 

Remember that to one-loop order the expression for the form factor reads [26] 



FvM = 1 + g^knr - *mi)JM + "-f (V> + , (7) 

where the function J is defined in [26J and where f n = 92.2 MeV denotes the 
pion decay constant. The small kaon loop contribution is taken care of by the 
use of the empirical value of the pion squared radius. 

The full expression for the form factor, Eq. fl6]), and its one-loop counterpart, 
Eq. (J7J), which enters the 1-loop ChPT prediction for 77 W 7r + 7r~7 of Ref. [3] 1 l 
are compared with data in Fig.[TJ While by construction both curves corre- 
spond to the same pion radius, the two descriptions start to deviate already 



1 Actually, Ref. [3] discusses the radiative two-pion decays of the octet and singlet 
states, 77s and r\\. 
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Fig. 2. Experimental data and error weighted fits for 77 (left, data are from Ref. [27] 
(filled squares) and Ref. [25] (open circles)) and rf (right, data are from Ref. [29] ) 
to 7r + 7r~7 according to Eqs. (UJ) and with s w = m (,) (rn (,) — 2Ey). The thin 
(green) line in the right panel denotes the possible impact of p — u mixing under 
the assumption that it appears here with the same strength as in in- 
visibly at values of s n7T as low as 0.09 GeV 2 . The kinematic range covered 
in the radiative rj and rj' decays spans from 0.077 GeV 2 to 0.301 GeV 2 and 
0.918 GeV 2 , respectively. The two ranges are indicated in the figure by the 
thick bars. While one might hope to be able to describe the form factor in the 
kinematic range for the 77 within ChPT at sufficiently high orders (e.g. below 
0.25 GeV 2 the 2-loop order seems to be sufficient [T^|T5] ). clearly a description 
for the rj' is impossible with any perturbative series. 

As we will demonstrate below: once the non-perturbative part of the transition 
amplitude in form of the pion vector form factor is divided out, what remains 
can be described by a polynomial in s n7T of low order that can be analyzed 
within U (3) extended ChPT, even in the case of the radiative 77' decay. 

3. We now turn to the evaluation of the full decay amplitudes. As mentioned 
above the first step is to analyze the data for both the total decay rates as 
well as the spectra with a fit of the parameters A and a of Eq. ([TJ). Note that 
the experiments discussed below confirmed that the pion pair is in the vector 
isovector channel. The resulting fits to the spectra of the WASA-at-COSY 
collaboration [27] for the 77 case and of the CRYSTAL BARREL collaboration 
[2U| for the 77' case are shown in Fig. [2J left and right panel, respectively. From 
the error weighted fits, we extract values of 

a = (1.96 ± 0.27 ± 0.02) GeV" 2 ; a' = (1.80 ± 0.49 ± 0.04) GeV -2 (8) 

where the parameter extracted from the data on the 77' appears as primed. The 
first and second uncertainty originate from the fit to the data on 77^ — > 7r + n~j 
and from that of the pion vector form factor, respectively. The a parameter 
was also determined directly by the WASA-at-COSY collaboration quoting in 
addition a systematic uncertainty of 0.59 GeV" 2 . Since the CRYSTAL BAR- 
REL data points include systematic uncertainties, the uncertainty of the a' 
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Fig. 3. Plot of P(s nn ), defined in Eqs. ([I]) and ([2]), as extracted from data for 
r\ 7r + 7r _ 7 [27] (left panel) and rj — > 7r + 7r~7 [29J (right panel). 

value should include both statistical and systematic uncertainties. We also 
studied other data sets for 77 and 77'. Concerning the former decay, Gorm- 
ley et al. [28] provides a = (1.8 ± 0.4) GeV -2 while Layter et al. [30] gives 
a = (—0.9 ± 0.1) GeV -2 . The acceptance correction of these old experiments 
was derived from the specified dT /dE 1 distributions, respectively, under the 
assumption that the pertinent matrix element is the simplest gauge invariant 
one (corresponding here to P (s W7T ) and F v (s nn ) equal to one). The Layter et 
al. result seems to be inconsistent both with WASA [27] and Gormley et al. 
[28] . However, from the information provided in those old experimental pa- 
pers it is impossible to evaluate systematic uncertainties. In case of the r/', we 
obtain a' = (2.7 ± 1.0) GeV" 2 from the data of the GAMS-200 collaboration 
[3Tj . which is larger, but within error bars consistent with the value listed 
above. Hence, in the following, we use the values given in Eq. (jHJ). 

In the pion vector form factor p — u mixing shows up as a quite spectacular 
effect. In case of the 7/ decay, as a result of the behavior of the rate at 
low values of E 7 , the effect is a lot smaller: if we take the mixing effect with 
the same strength as it appears in F v using the prescription of Ref. [18], the 
impact of the mixing on the r( spectrum is very moderate — see the wiggly 
(green) line in the left panel of Fig. [2J A fit to the rj' spectrum including the 
mixing as shown shifts the value of a' upwards by 0.3 GeV -2 , well within 
errors. 

Instead of looking at the data themselves it is illustrative to extract from data 
directly the polynomials P(s n7T ). These are shown for both radiative r/ and rj' 
decays in the left and right panel of Fig. [3], respectively. Here one clearly sees 
that the residual s n7T dependence for both transition amplitudes — once the 
pion form factor and the phase space are divided out — has a linear behavior to 
a very good approximation. The statement is further corroborated by the fact 
that any additional quadratic term to the linear polynomial with coefficients 
as specified in Eq. ([8]) is compatible with zero: = (0.21 ±0.67 ±0.06) GeV -4 
and (3' = (0.04 ± 0.36 ± 0.03) GeV -4 . This appears reassuring, although it 
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came as a surprise that even for the 77' a first-order polynomial is sufficient. 
The origin of this might be in the current quality of the data which is best in 
the region of large values of E 1 which corresponds to moderate values of s n7T 
- this is the region where the chiral expansion is expected to converge (once 
resonance effects are taken out). This can also be seen in Fig.[3l right panel: 
clearly the fit is dominated by values of < 0.6 GeV 2 (this corresponds to 
pion relative momenta of at most 360 MeV) that are still reasonably smaller 
than the typical hadronic scale of order of 1 GeV, which sets the expected 
range of convergence for the chiral expansion once the vector pion form factor 
is taken out. 

The normalization of the data used above cannot be employed to fix the 
prefactor for they are given in arbitrary units only. However, once the 
slope parameters are fixed we can use the experimentally measured partial 
widths [32], = (59.8 ± 3.8) eV and = (57.0 ± 2.8) keV, to 

extract 

5 = -0.22 ±0.04 ; 5' = -0.40 ± 0.09 , (9) 

where we used the definition = A^\l + 5^), with Aq for the transition 
strength in the chiral limitEi which is fixed for the transitions at hand by the 
chiral (box) anomaly [3117] . 

It should be stressed that the values extracted for are a lot more reliable 
at this stage than those extracted for \ since the former are differential 
quantities while the latter are integrated ones, being quite sensitive to the 
shape of the spectrum also in the regime where we do not have high quality 
data yet. Thus, one might expect that, once these become available also for the 
rj decay spectra at large values of s n7T , a linear polynomial is insufficient for 
P(s n7T ). This should not change \ however, it might significantly influence 
the integrated rate and therefore the values of 

4. In order to perform the matching of the decay amplitude to one-loop ChPT 
we now replace Fy in Eq. ([TJ by its one-loop expression, Eq. ((Tj), and expand 
the result to first order in s n7T . The resulting expression can then be equated 
with the corresponding one from one-loop ChPT, see appendix |A] This pro- 
cedure gives a relation between the phenomenological parameters and 5^ 
and the U(3) ChPT low-energy parameters a[ , a^ K \ a±~ n \ and 02- 



2 Here and in the following we will use the term 'chiral limit' in a somewhat loose 
sense, since we insert in Aq the physical value of the pion decay constant f w and 
assume that /g 7^ f n for the octet pseudoscalar decay constant, etc. 
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Especially one finds from the independent terms for the 77-decay 
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and from the s nn dependent terms for the 77-decay 
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Here denotes the kaon mass. In line with previous investigations we iden- 
tify the scale \x with m p . Furthermore we introduced the abbreviations 



Af, 
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4v^vr2/3 fl 



with e the unit of electric charge, a = cos(0) and b = — sin(6*) in terms of 
9 ~ —20° as the value of the 77 — 77' mixing angle and /§ and /1 the values 
of the octet and singlet pseudoscalar decay constants - all three values as 
specified as in Refs. [6fT] , see also Ref . [3] . Since our aim here is the comparison 
with the existing radiative decay results of the above mentioned references, we 
still follow the old parameterization in terms of two (octet and singlet) decay 
constants and one 77 — 77' mixing angle. This is in contrast to more modern 
parameterizations in terms of either two octet-singlet decay constants and two 
mixing angles, which follow from the defining matrix elements of the octet and 
singlet axial- vector current [33], or in terms of strange and non-strange decay 
constants and only one mixing angle, see e.g. Refs. [34|35f36] . 

The transition strength in the chiral limit is given by A = A 8 +Ai. To arrive 
at the corresponding expressions for the 77' decays, a, 5, A , a and b need to be 
replaced by their primed counter parts, especially a' = sin(0) and b' = cos(6 l ). 
In Eq. pip we have only kept the leading term of the kaon loop contributions, 
which turns out to represent the complete expression to high accuracy. 

In the case of the independent terms, we observe that two parameters 
from the phenomenological expression, 8 and 5', have to be matched onto four 
parameters from U(3) extended ChPT, af n \ af K \ and a[ . Thus the 
matching does not provide any constraint. The situation is different for the 
s nn dependent terms, as the two parameters a and a' have to be matched 
to the single parameter 02- To justify this statement we needed to ignore 
s nn dependent counter terms subleading in N c , the number of colors. Note, 
however, that these neglected terms are used to cancel some of the additional 



9 



divergences in the one-loop approximation of the rjx — > tt + ti 7 decay [3]. 
Extracting ci2 from the rj and the 7/ decays gives 

a\ = 9.70 ± 0.7 ; 4 = 9.23 ± 1.4 , 

respectively. Thus we indeed find that the data sets for the rj and 77' decays 
are consistent with the assumption that the only non-perturbative part of the 
amplitude originates from the pion vector form factor. 

5. In the following, by invoking chiral Ward identity and large N c arguments, 
we will give a physical interpretation to essentially half of the parameter \ 
When the hidden gauge approach [3"Tf3"8] to low-energy hadron physics was 
first applied to the anomalous sector, it was already found out that the so- 
called complete VMD (for both the triangle and the box anomaly sector) 
was unsustainable in this approach [39], see also the reviews |40f41f42] : while 
the non-anomalous and the triangle-anomaly sectors were fully compatible 
with VMD, the description of the box-anomaly-induced decays gave satisfying 
predictions (see e.g. the u — > nun decay) only if the VMD triangle anomaly 
terms were supported by a point-vertex (contact term) involving a photon 
and three pseudoscalars. Cohen [43] was the first (see also [44]) to point out 
that the chiral Ward identity implies that both the chiral triangle and the 
box anomaly contribute to the ^nnn processes and the r/ ) — > nn^ decays 
away from the chiral limit. With the help of the low-energy theorem for the 
7 — > tttttt amplitude of Ref . [45] , he argued that the total amplitude for these 
(in the chiral limit by the box-anomaly induced) processes can be decomposed 
as 

^tot _ 3^VVA _ l^VAAA 

where the triangle or VVA-type amplitude [V: vector, A: axialvector] , under the 
assumption that one of the vectors of the VVA-type amplitude subsequently 
decays into two pions, and the box or VAAA-type amplitude contribute with 
a relative weig ht of | : -§. Still there is an ongoing debate whether these 
VAAA-type contact terms are empirically needed or not, see e.g. [46|5] . 

The vector pion form factor Fv(s n7T ), which - as mentioned before - contains 
both the Born terms and final state nn interactions, can be represented by 
one of the V legs in the triangle anomaly VVA vertex (the other V leg denotes 
the outgoing (isovector) photon, whereas the A leg stands for the decaying 
pseudoscalar 77 under PCAC). In other words the VVA vertex 

2*A rj — >7r7T7(0)-PV(s 7rYr )-f > (s 7r7r ) (13) 

(with P(s n7T ) = 1 + as n7T + 0(sl n ) a polynomial, similar to P(s 7r7r ) in Eq. ([1])) 
contains both tree-level contributions of the vector form factor, which are 
of leading order in the 1/N C expansion, and genuine loop contributions of 
Fy, which are suppressed by at least one factor of 1/N C and are therefore 
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subleading_i|. In the case of the VAAA box anomaly, however, the two remaining 
AA legs stand for two separate pions which either are already the final pions 
or which still rescatter. The first case, which is the leading- N c contribution of 
the VAAA process (times the 0{N^ 1 ^ 2 ) scaling of the prefactor ^^^(O)), 
would correspond in the VVA scenario just to the replacement Fv(s nn ) — > 1, 
i.e. to the trivial term in the Taylor-expansion of the vector form factor ([6]), 
while the second (rescattering) case, since it necessarily involves an additional 
four-pion vertex, is represented by the subleading terms in the 1/N C expansion 
of Fy(s nn ). The leading-order terms of the VVA- and VAAA-type amplitudes 
have to have the same N c scaling, since both are tree-level processes and since 
the initial state (r/-') and the final state (7r + 7r~7) are the same, respectively. 
However, only the former can run through a vector meson pole at tree-level. 
The total VAAA result is therefore 

- |A->-7(0)e-^ (r2>s ^Fy(^)P(^ 7r ) . (14) 

Here the coefficient (r 2 ) is the leading N® contribution of the mean square 
charge radius of the pion, which to 0(m 2 ) in ChPT is given by (r 2 ) = (/6 — 
l)/(167r 2 /2) ; S ee e.g. Ref. [20], such that (r 2 ) = (k-l" S ) I '(16vr 2 / 2 ). Remember 
that Iq itself is of order N c and that / w scales as \/W c . The exponential term 
exp(— |(r 2 )s 7r7r ) in the VAAA vertex therefore serves to remove the leading N® 
contributions in Fy(s nn ), such that - with exception of the trivial term - only 
the subleading ones are left over (multiplying — ^A^^^O)). Summing the 
VAA and VAAA induced contributions to the rf^ — > nn^ decay amplitude, 
i.e. ({TBI) and (fl4"|) . and Taylor-expanding the exponential factor of (fl4"|) and 
P(s mr ), we therefore get the total amplitude 



7r7ry(^7r7r 



■77— >7T7Vy 



(0)F, 



V{S? 



+ (^(r 2 ) + a) s n7T + 0(s 



2 y 

7T7VJ 



(15) 



Thus, from this point of view, the coefficient a comprises of two terms: the 
leading N® term of the subtraction constant of the vector form factor, | (r 2 )/6 ~ 
/6/(1927r 2 / 2 ) Rs 1.03 GeV^ 2 [20], and the remainder a. Comparing to the num- 
bers of Eq. (JSJ), we find in this case that about half of the value of can be 
interpreted as coming from (r 2 ). 

6. Finally, we will link our method to earlier studies and compare results with 
vector meson dominance considerations in general. Inspired by an N/D anal- 
ysis of the process 7 — > ixixtx (see also [8ll4"6] ). Holstein and Venugopal [6]T7] 



3 Note that the prefactor ^l j; _ il7r7r7 (0) contains the coefficients N c / f% ~ 0(N C 1//2 ) 
and therefore determines the leading overall A^ c scaling behavior of the total 7/ > — > 
7T7T7 decay amplitude. It is of course the same for the VVA and the VAAA case. 
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constructed a precursor to the relation ([T]). They started from an ansatz con- 
taining a contact term as well as a rescattering term parameterizing the tttt 
final state interactions, which reads in the notation used above 

A Tl ^ n7Ty (s nn ) = A^nir-y [1 - c + c(l + as nn )F v (s nn )] . (16) 

Here c and a are free real-valued parameters and .A^ is the amplitude in 
the chiral limit. By matching this ansatz - both at order 0(p 6 ) - to one- loop 
chiral perturbation theory [3J (more precisely, to the coefficient in Ref. [3J of 
the standard one- loop function J(s n7T ), cf. Eq. ([7J) on the one hand and to 
VMD (as in Ref. [4]) on the other hand, these free parameters were fixed in 
Refs. [6f7] to c = 1 via ChPT and to a = l/(2mV) via VMD, respectively. In 
light of the discussion presented above and as implicitly stated in Refs. [37], 
c = 1 {i.e. no contact term) appears as a necessity from dispersion theory. 
Note that Ko and Truong [5] derived the same expression from unitarity and 
the above-discussed Ward identity constraints. If the parameter a in Refs. [6f7] 
were not extracted from VMD or a 1 = were not implicitly assumed as in 
Ref. [8] and the amplitude ^^^(O) were not fixed to its chiral limit value 
A®^ nnry , then the resulting relations of Refs. [617] and of Ref. [5] would have 
had the form of Eq. flTJ, for P(s n7T ) expanded to first order. 

Let us now compare to vector meson dominance in general: in its most simple 
form one would just get a = 0, at variance with data, see (jSJ). If, following 
Ref. [43], the chiral Ward identities are implemented, then the simplest sce- 
nario would correspond to a = instead of a = 0, a result at the edge of 
being consistent with current data. This is what was used in Refs. [6H7II8] . 
Finally note that in all evaluations that modify the r\ — > 7r + 7r~7 decay ampli- 
tude of the chiral limit solely by s Tn dependent terms, as is the case in vector 
meson dominance models, it is impossible to simultaneously predict the ex- 
perimentally measured shape of the distribution and the empirical branching 
ratio. 

7. To summarize, the rf^ — > 7r + 7r~7 decay amplitude, here assumed to be 
p-wave dominated, factorizes into a universal part and a reaction specific part 
after the trivial (point-particle) kinematics is removed. If higher partial waves 
cannot be neglected, their contribution should be added suitably to the above- 
mentioned p-wave amplitude. The universal part, which applies to all p-wave 
dominated radiative decays with an isovector 7t + 7t~ pion pair in the final state, 
is given by the well established vector pion form factor Fy. The reaction spe- 
cific part can be parametrized as an expansion in the invariant mass s n7T of 
the pion pair: P (s nn ) times the normalization factor A — as long as left-hand 
cut contributions are suppressed, as is the case for the reactions under con- 
sideration. The expansion can then be treated perturbatively, since unitarity 
and analyticity of the final-state interaction dictate that the pion form fac- 
tor already takes care of the non-perturbative aspects of the tttt unitary cut. 
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Moreover, the perturbative expansion allows a systematic comparison with 
ChPT predictions. Especially, the case with ?/ decays shows a way to extend 
methodically perturbative calculations also to the resonance region. In this 
sense the presented approach is indeed model-independent. 

For the description of the present world data it is sufficient to consider only 
the linear term (the ) parameter) in the expansion ofP (s nw ) = 1 + ^s n7T + 
0(s^ n ) — both for the r\ and rj decays. The extracted values of this parameter 
for the existing experiments are a = (1.96 ± 0.27 sta t ± l-00 syst ) GeV~ 2 and 
a' = (1.80 ± 0.49 tot ) GeV -2 . The data of Gormley et al. %ZE\ and the new 
data of the WASA-at-COSY collaboration [27] are consistent (both support 
an a parameter of the order 2 GeV -2 ), whereas the data of Layter et al. [3U] 
indicate an a value below zero. However, the old experiments (Gormley and 
Layter) do not provide estimates of the systematic uncertainties. 

We propose to use the introduced a and a' coefficients or, more generally, 
the polynomial P(s nn ) to parametrize and compare - including the pertinent 
statistical and systematical uncertainties - the spectra of all past and future 
radiative decay experiments with only one (iso)vector p-wave 7r + 7r - pair in 
the final state. New experimental data on the 77 — > 7r + 7r~7 decay distributions 
from KLOE are at the final stages of the analysis jUJ and will be released soon. 
Meanwhile also WASA-at-COSY, CLAS and BES-III have collected new large 
data sets of respectively 77 and 7/ decays that will enable us to check if further 
terms in the expansion of P(s 7r7r ) are necessary. 
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A Explicit expressions from one— loop U(3) extended ChPT 

The ChPT amplitude employed in the matching (cf. Eqs. (flO| [TT]) ) has been 
obtained as follows. The divergences occurring in the regularized one-loop 
expression [3] are absorbed by anomalous 0(p 6 ) counter terms - a complete 
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list of which is presented in [4*8f4"9] in the SU(3) case. We follow Ref. [9] where 
also some U(3) extensions of these terms can be found. However, we neglect the 
pure singlet terms containing the coefficients W13 and Wn of Ref. [9], which 
posses an s^-dependence, but have one more trace than their SU(3) counter 
parts, such that they are only subleading in the 1/N C expansion. Furthermore, 
we use the Gell-Mann-Okubo formulae to dispose the rj [50] and 77' [51] masses, 
which is correct to 0(p 6 ). In summary, we get the ChPT amplitude 



A 

with 



ChPT 



A 8 ■ C va + At ■ Cje^p {e*Y{pJ v { P+ Y{p-Y (A.l) 



C m = l + C l r ps + 



1 



'As 



(8tt) 2 , (8K) 2 , 

a\ 'm n + a\ 'm K + a 2 s 7 



a 



m 



/^yloops 1 



32tt 2 /. 2 

1 r (Itt) 2 , (IK") 2 , 

a\ 'm n + a\ 'm K + a 2 s 



(A.2) 
(A.3) 



Here the coefficients C^° ops comprehend the finite loop contributions (see [3]), 
while the remaining constants comprise the low-energy constants (LECs) of 
Ref. [9] in the following way: 



,(8*0 



J 1K ) 



-Aw 



2 -(0) 



K 



to>P + |tDg> 



4^ 0) + 



2^ 

„T,(°) 



Z 7 T;(°) 



124 0) 



5*;}? 



2«$ 



6u>i + 4w 



t,(°) 
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a 2 = k 



2^? + «>g> 



(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 



with k = 2 11 7r 4 /^. Hence the momentum dependent terms of t] 8 and rjx involve 
only one linear combination of the LECs (cf. a 2 ), whereas the mass terms are 
governed by four linear combinations. 



Note that the structure of C m and C Vl displayed in Eqs. (1A.2|) and (1A.3[) 
can also be derived from [4"8H4"9] . if the Gell-Mann-Okubo formulae and the 
pertinent U(3) extension are applied as above. 
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